Electric-field controlled spin in bilayer triangular graphene quantum dots 
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We present theoretical results based on mean-field and exact many-body approaches showing that 
in bilayer triangular graphene quantum dots with zigzag edges the magnetism can be controlled by 
an external vertical electric- field. We demonstrate that without electric field the spins of the two 
layers of the quantum dot interact ferromagnetically. At a critical value of the electric- field, the 
total spin of the bilayer structure can be turned off or reduced to a single localized spin, a qubit 
isolated from contacts and free from interaction with nuclear spins. 



Graphene exhibits unusual electronic properties [l|- 
including relativist ic nature of quasi-particles, sublat- 
tice pseudospin, and zero energy bandgap. The appli- 
cation of graphene for logic devices requires opening of 
the bandgapwhich can be achieved by either size quan- 
tization [Thllq. or chemical modification [19], or bring- 
ing in a second layer and applying an external electric- 
field [2Q|-r28.] , or size quantization in bilayer graphene 
nanostructures [29l - [32t . In particular, when graphene is 
reduced to a triangular quantum dot with zigzag edges [sj, 
sublattice symmetry is broken, size dependent energy gap 
is open, and a band of degenerate states at the Fermi level 
is created leading to finite macroscopic spin polarization 
[m-TlSj. This allows simultaneous size, shape and edge 
engineering of magnetic, electrical, and optical properties 
in a single material-graphene[14j. 

In this work, we investigate electronic and magnetic 
properties of bilayer triangular graphene quantum dots 
with zigzag edges under external vertical electric-field. 
We show that the magnetic moment of bilayer triangular 
graphene quantum dots can be controlled by the ver- 
tical electric-field. Without the electric-field, the mag- 
netic moments of the two layers are shown to be coupled 
ferromagnetically. Using configuration interaction and 
mean-field calculations based on tight-binding model, 
we demonstrate that the ferromagnetism can be either 
turned off or reduced to a single electron/hole spin. The 
single electron spin is hence isolated in a charge neutral 
structure by the application of an electric-field, indepen- 
dent of the size of the quantum dot and without deco- 
herence due to contacts. The electric-field control of the 
ferromagnetism [33] and isolation of a single spin opens 
new applications in spintronics and quantum information 
processing jsl - lssj . 

Figure la shows two possibilities for building a bilayer 
triangular graphene quantum dot (BQD) using two sin- 
gle layer triangular quantum dots (TGQD) of comparable 
sizes, with zigzag edges. We consider AB Bernal stack- 
ing, where the A sublattice of the top layer (A2, shown 
in blue color) is on top of the B sublattice of the bot- 
tom layer (Bl, shown in red). On the left hand side, the 
two TGQDs are of the same size. In this configuration. 
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FIG. 1: (Color online) Bilayer triangular graphene quantum 
dot with zigzag edges, (a) structure constructed using two 
single layer quantum dots with equal sizes (on the left) and 
different sizes (on the right), (b) Isosurface plot of the spin 
density — Pi for with and without electric- field E, obtained 
from configuration interaction calculations. 



however, not all the A2 atoms have a Bl partner as re- 
quired by Bernal stacking. A more natural configuration 
choice is shown on the right hand side of Fig. la. The 
top layer triangle has its fioating atoms removed, mak- 
ing it smaller than the bottom layer triangle. Such a 
bilayer construction has the interesting property of hav- 
ing an odd number of degenerate states at the Fermi 
level (zero-energy shell of states) independent of its size, 
allowing to isolate a single spin in a charge neutral struc- 
ture and hence isolated from the contacts. This is il- 
lustrated in Fig. lb, where the spin density isosurfaces 
are shown for zero electric-field (left hand side) and fi- 
nite electric-field (right hand side), as obtained from our 
tight-binding based configuration interaction calculations 
explained in detail below. When the electric-field is off, 
both layers have a finite magnetic moment, as in single 
layer triangles [12l-[l8|. differing by one spin due to the size 
difference of the two triangles. The magnetic moments 
of the two layers are coupled ferromagnetically. When 
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FIG. 2: (Color online) Single particle tight-binding spectrum, 
(a) The bilayer quantum dot consisting of 1195 atoms has 43 
zero-energy states, (b) When an electric-field is applied, the 
degeneracy between the 21 top layer zero-energy states and 
22 bottom layer zero-energy states is lifted. 



a sufficiently high electric field is applied, electrons from 
the lower layer reduce their energy by flipping their spin 
and transfering to the top layer filling up all the available 
spin up and down zero-energy states, leaving behind one 
single spin. 

We have confirmed the above picture through calcula- 
tions with varying levels of accuracy. As a first step, we 
diagonalize the tight-binding Hamiltonian given by 



(1) 



where the operator cj^ creates an electron on a Pz orbital 
on site i with spin a. The tight-binding parameters r^j 
are fixed to their bulk values t = —2.8 eV for in-plane 
nearest neighbours i and j and t± = 0.4 eV for inter-layer 
hopping. The effect of the potential difference induced 
by an external perpendicular electric-field E is taken into 
account through Vi = —AV/2 for the bottom layer atoms 
and Vi = AV/2 for the top layer atoms. Figure 2a shows 
the energy spectrum near the Fermi level for AV = 
for a BQD consisting of 622 atoms in the bottom and 
573 atoms in the top layer. If we take t± = 0, the two 
triangles are decoupled and we find 22 zero-energy states 
in the bottom layer and 21 zero-energy states in the top 
layer, for a total of 43 zero-energy states^onsistent with 



previous work on single layer TGQDs|12l-ll8j. Turning 



on t± to 0.4 eV does not affect the zero-energy states. 
The effect of applying an electric field, e.g. AV = 0.4 
eV, is shown in Fig. 2b. The energy of the 21 zero-energy 
states corresponding to the top layer is pushed up by 0.4 
eV with respect to the bottom layer zero-energy states. 
Note that the bottom layer zero-energy states do not ex- 
perience any dispersion unlike the top layer zero-energy 
states. This is due to the fact that they lie strictly on Al 
sites which are not coupled to the top layer, whereas the 
top layer zero-energy states, which lie on B2, do couple 
to the bottom layer. The ability of controlling the rel- 
ative position of zero-energy states gives an interesting 
opportunity to control the charge and spin of the zero- 
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FIG. 3: (Color online) Meanfield Hubbard results. Energies 
of lowest energy states with different total spin S as a function 
of potential difference AVc between the layers, with respect 
to the ferromagnetic configuration Smax = 4.5 for a system 
of 107 atoms and 9 zero-energy states. 



energy states, expected to be spin polarized for AV = 
according to Lieb's theorem[39]. 

The magnetic properties of BQD structures as a func- 
tion of applied electric field can be studied by solving 
the Hubbard model within the self-consistent mean-field 
approach: 



H 



MF 



^ ViUi, + U - -)n,,(2) 



where U is the on-site Hubbard term taken to be 2.75 
eV here. First, we study a small BQD of 107 atoms 
and 9 zero-energy states. Figure 3 shows the energies 
for different total spin S with respect to the energy of 
the ferromagnetic configuration, = 9/2. At AV = 0, 
the degenerate band of zero-energy states is polarized: 
all 9 electrons occupying the 9 zero-energy states have 
their spins aligned ferromagnetically. Although the first 
excited state obtained from the Hubbard model is an- 
tiferromagnetic with the total spin of the bottom layer 
opposite to the total spin of the top layer, a full treat- 
ment of the correlation effects shows that (see below) low 
lying excited states have more complex spin structures. 
The Hubbard model is, however, useful for estimating the 
critical value Vc where phase transition occurs. As AV is 
increased, the electrons lying on the bottom layer zero- 
energy states are forced to flip their spin and tunnel to 
the top layer zero-energy states. At around AVc = 0-55 
eV such charge transfers occur abruptly, leading to a de- 
crease of the total spin of the system. As a result, all top 
layer zero-energy states become doubly occupied, leaving 
exactly one single spin in the bottom layer zero-energy 
states. We note that one can also isolate a single hole 
spin in the bottom layer by applying a reverse electric 
field, thus pushing the electrons from the top layer to 
the bottom layer, occupying all states except one. It is 
thus possible to isolate and manipulate a single electron 
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FIG. 4: (Color online) Size dependence of ferromagnetic- 
antiferromagnetic transition, (a) FM-AFM energy difference 
E{Smin) — E{Smax) as a function of potential difference AV 
between layers, up to A/" = 1507 atoms, (b) For AV = 0, 
the FM-AFM energy gap per number of side atoms Ngide ap- 
proaches 14.3 meV (top panel). Critical value AVc where the 
transition occurs as a function of number of atoms N (bottom 
panel) . 



or hole spin in a neutral BQD by applying an external 
electric field. 

The procedure of isolating single electron or hole 
should occur regardless of the size of the system since the 
top layer has always one less zero-energy state than the 
bottom one. In order to investigate the size dependence, 
in Fig. 4a we show the energy difference between the ferro- 
magnetic and antiferromagnetic (FM-AFM) states calcu- 
lated in the mean-field Hubbard approximation as func- 
tion of applied voltage for several sizes up to 1507 atoms. 
We note that, due to the unusually high degeneracy of the 
states, self-consistent iterations occasionally get trapped 
in a local energy minima. We have thus repeated the cal- 
culations several times using different initial conditions 
and/or convergence schemes to assure that the correct 
ground state was reached. As expected, at AV = 0, the 
FM-AFM gap increases with the size of the system N. 
In fact, the FM-AFM gap energy per Ngide^ the number 
of side atoms on the top layer (Ngide) (equal the number 
of inter-layer bonds on the edges), approaches a constant 
value of 14.3 meV as shown in Fig. 4b. However, the FM- 
AF transition voltage Vc decreases with the system size 
as can be seen from Fig. 4c. For the largest system size 
studied, N = 1507, we obtain AVc = 0.345 eV, which 
corresponds to an electrical field of ~ 1 V/nm, a value 
within experimental range [24]. 

We now test the predictions of the mean-field Hubbard 
model by including long range interaction and correla- 
tion effects, which is computationally feasible for the sys- 
tem size studied in Fig. 3. We first solve a Hartree-Fock 
generalization of the Eq.2, but with empty zero-energy 
states [13|]: 



FIG. 5: (Color online) Configuration interaction results, (a) 
Energies of states for different total spin 5* as a function of 
potential difference AV between the layers, with respect to 
the ferromagnetic configuration Smax =4.5 for the same sys- 
tem as in Fig. 3. (b) Ground state total spin 5* as a function 
of Ay for the charge neutral system, — e charged system, and 
+e charged system. 
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where the tight-binding term now includes the intra- 
layer next-nearest neighbour hopping t = —0.1 eV, inter- 
layer next-nearest hoppings 73 = 0.3 eV and 74 = 0.15 
eV[27]. The terms p and p^'^^^ are quantum dot and 
bulk density matrices, respectively [13|. The two-body 
Coulomb matrix elements {ij\Vee\kl) computed using 
Slater pz orbit als include on-site interactions, all scatter- 
ing and exchange terms within next-nearest neighbors, 
and all long range direct interactions. We have pre- 
viously tested the validity of our approach by compar- 
ing to density functional calculations and obtained good 
agreement [13]. After diagonalizing the Hartree-Fock 
Hamiltonian, we obtain Hartree-Fock quasi-particles de- 
noted by the creation operator with eigenvalues 
and eigenfunctions \p). We can now fill the new quasi- 
particle zero-energy states with electrons and solve the 
many-body Hamiltonian given by 

pa pqrsaa' 

Figure 5a shows the evolution of spin states for the 
same system studied in Fig. 3, but now obtained by diag- 
onalizing the many-body Hamiltonian. We observe two 
main differences from the mean-field Hubbard results: (i) 
At low AV, the antiferromagnetic configuration S = 1/2 
is no longer the first excited state. Although the AFM- 
FM energy gap is still comparable to the Hubbard result, 
other spin excitations are now closer to the ground state 
due to correlation effects, (ii) The ground state spin tran- 
sitions do not occur as abruptly as in a Hubbard model. 



4 



As shown in Fig. 5b, the total spin of the charge neutral 
system (red color online, solid line with circles) evolves 
towards the minimum spin state S = 1/2 gradually be- 
tween /SV = 0.5 — 1.0 eV. This is mainly due to long 
range interactions. As the electrons are transferred one 
by one from the top layer into the bottom layer, they 
leave a positively charged hole behind which makes it 
harder to transfer more electrons. We note that this be- 
havior of gradual spin transition is also obtained within a 
mean-field Hubbard model with long range interactions 
included (not shown). Finally, in Fig. 5b, we also study 
the effect of charging the BQD system. The asymmetry 
between magnetic moment of -he and — e charged sys- 
tems refiects the correlation induced spin depolarization 
process that occurs in single layer TGQDs as discussed 
in our previous works[13j. 

In summary, we demonstrate that the graphene bilayer 
triangular quantum dots exhibit a shell of degenerate 
states at the Fermi level. At half filling, the shell is max- 
imally spin polarized. By the application of a vertical 
electric field the total spin of the bilayer structure can be 
turned off or reduced to a single localized spin, a qubit 
isolated from contacts and free from interaction with nu- 
clear spins. This opens new possibilities in the area of 
spintronics and quantum information processing. 
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